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In this work, we have considered FRW model of the universe for Brans-Dicke (BD) 
theory with BD scalar field as a Chameleon field. First we have transformed the field 
equations and conservation equation from Jordan's frame to Einstein's frame. We 
have shown in presence of variable modified Chaplygin gas, the potential function 
V and another analytic function / always increase with respect to BD-Chameleon 
scalar field (p. 
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1. INTRODUCTION 

The present accelerated expansion of the universe are strongly believed under some 
observations of the luminosity-redshift relation for type-la supernovae [1, 2], measurements 
of CMBR, power spectrum of mass perturbation etc. The responsible candidate for the 
violation of strong energy condition is known as dark energy. Another simplest alternative 
candidate which includes the scalar field in addition to the tensor field in general relativity is 
Brans-Dicke (BD) theory. Recently, enormous interest in BD theory [3] can be found to the 
discovery by La et al [4] that the use of BD theory in place of general relativity can upgrade 
the problem of inflationary cosmology. This is possible as it arises naturally as the low energy 
hmit due to the interaction of the BD scalar field. The field equations in this case are modified 
due to the coupling between matter Lagrangian and scalar field. In BD theory, there exists 
one dimensionless parameter u and the effective gravitational constant G which is inversely 
proportional to the scalar field (j). Also, BD theory or its modifications become an important 
theory to explain the present cosmic acceleration [5] with some correction of Newtonian 
weak-field limit [6]. Einstein gravity can be reconstructed from it considering the condition 
a; — > oo [7]. Several authors have done extensive amount of work on BD theory and gen- 
eralized BD theory to solve different cosmological problems of inflation, quintessence etc [8 - 12]. 

There are different dark energy models which are basically considered as a scalar field, 
suggest a large correction to the Newton's law. These are nearly massless scalar field 
and if exists on earth, should have been detected in local test of equivalence principle 
and as a fifth force. Recently, a different approach is introduced in general relativity by 
Khoury and Weltman [13], where a non-minimal interaction occurs between the quintessence 
scalar field and matter sector rather than the geometry and this interaction is introduced 
through an interference term in the action known as Chameleon scalar field. The important 
feature of this field is that due to the coupling to matter it obtains an effective potential 
i.e., an effective mass depends quite extremely on the background energy density. Many au- 
thors have been studied on Chameleon field in different aspects of accelerating universe [14 - 21]. 

A well known candidate for Q-matter is the exotic type of fluid named Chaplygin gas 
obeying the EOS p = —— {B > Q) [22], where p and p are respectively the pressure and 
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energy density. A generalized version of this Q-matter having the form of EOS given by 
P = -^(0<q;<1) [23,24] and recently it was modified to the form p ^ Ap - ^ {A > 0) 
[25-27] , which is known as modified Chaplygin Gas. This model can represents the evaluation 
of the Universe starting from the radiation era to the ACDM model. There is one stage, when 
the pressure vanishes and the matter content is equivalent to a pure dust. Recently another 
modification have been made by Guo and Jhang [28]. They proposed variable Chaplygin gas 
model (VCG) where S is a -|-ve function of the cosmological scale factor a. This B{a) is 
related to the scalar potential if we take the Chaplygin gas Born-Infeld scalar field [29] . 

From these three models (BD theory. Chameleon field and VCG) we have construct our 
paper as follows: BD theory or its modifications have recently described as a driver of the 
present cosmic acceleration. It had been shown in [5] that BD theory can potentially generate 
sufficient acceleration in the matter dominated era even without including the exotic Q-field 
[5]. But this has problems to estabhsh the transition from a decelerated to an accelerated 
phase. On the other hand, this transition is possible when a Chameleon scalar field have a 
non-minimal coupling with dark matter [20,21]. Here we consider the BD scalar field as a 
Chameleon type scalar field and the matter Lagrangian is coupled with an arbitrary analytical 
function of the field. Initially, we write the field equation in most general scalar tensor theory 
of gravity but the field equations become complicated for the gravitational and scalar fields. 
So we transfer these equation in Einstein frame using conformal transformation [30,31], in 
which the scalar and tensor degrees of freedom do not mix. Now we consider the universe filled 
with variable modified Chaplygin gas (VMCG) [32] and in general scalar tensor gravity, we 
graphically analyze the potential of the BD-Chameleon scalar field with respect to the field. 
We also analyze the analytical function of the BD-Chameleon scalar field. 



II. BASIC EQUATIONS AND SOLUTIONS 



We consider here the Brans-Dicke scalar field (f) act as a Chameleon field and we consider 
the action for the chameleon scalar field in self- interacting Brans-Dicke (BD) theory [31] as: 
(choosing SttGo = c = 1) 
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where V{(f)) is the self-interacting potential for the BD-Chameleon scalar field 4> and cu is the 
BD parameter. Unlike the usual Einstein- Hilbert action, the matter lagrangian is considered 
as f{(f))Ljn-i where is some arbitrary analytical function of 0. This term makes pos- 

sible the non-minimal interaction between the cold dark matter and BD-Chameleon scalar field. 



The metric of a homogeneous and isotropic universe in the FRW model is 



dr^ 



+ r^((ie^ + sin2 Odc})'^ 



(2) 



1 — fcr^ 

where a{t) is the scale factor and A; (= 0, +1, — 1) is the curvature scalar. Now, in BD 
theory, the Einstein's field equations for the above space-time symmetry are 



rr2 k f{<f>) Vi(f)) 



6 02 
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(3) 



2- + H^ + — + + 2H^ + ^ 



a 2 02 ' " ' ^ ' 20 

and the wave equation for the BD-Chameleon scalar field is 



(4) 



(2a; + 3)0 



y(0)- 



1 dV" 
2'^d0. 



The energy conservation equation is given by 



'^z X 3, x/(0) 

a 4 /(0) 



Here we consider the variable modified Chaplygin gas with equation of state [32] 



(5) 



(6) 



p^Ap ,^>0,S>,0<q;<1 

pa 



(7) 



where B is not constant. In the BD theory, the effective gravitational constant is given by 
G — which is indeed not a constant. Now considering the the transformation 



(8) 
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so that the equations in the Jordan's frame reduce to the equations in the Einstein's frame. 
So let us define the new variables in Einstein's frame as in the following: 



ijj = In (J), a = (pa ,p = (pp, p = 0p, / = 0/, V = (t)V,dT = \j4>dt (9) 

Now the BD field equations (3), (4) and the conservation equation (6) can be reduced to the 
Einstein frame as 



a? ' a? V 12 



a" a'2 k /2cc; + 3\ ,,0 Vii^) 



and 



-p' + 3^(p + p) = -\{p - ?>PW - |(P + P) J (12) 

where, ' represents the derivative w.r.t. new time variable r. Also using the above trans- 
formation and considering the variable term B — Bq<P°'^^, the EOS of the variable modified 
Chaplygin gas (7) becomes 

P-AP-^ (13) 

pa 

where -Bo(> 0) is a constant. This EOS represents only the modified Chaplygin gas in 
the Einstein's frame. From the equations (12) and (13) we can construct a linear differential 
equation in z (where, p — 2;«+t) as 



_ + ,(„ + l)/3(l + ^)- + -(l + A)i. + — 



3^ ,A . X / ^ «' 3„ /' 3^ 
V' j = (a + 1) \3Bo- + -Bo J - -BotP' 



(14) 

For simplicity we take only the radiation case A = 1/3 and in this case the integrating 
factor (l.F.) becomes (a^/)""*"^. We also assume a^/ = 0", where n is a positive constant. After 
solving the above the differential equation, we get the solution of p as 



4 "I n(a + l) + l 



1 

a + l 



(15) 



where C is an integrating constant. Now let the transformed Hubble parameter is defined by 



^ — H then H' — aH^. Now we take a new assumption tp — ipoo^ where m is a positive 
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constant. Using this new assumption the solution of p becomes 

a + l 



-Boil- , 
4 n(a + l) + l 



(16) 



Now putting H"^ — Y and considering the previous all assumptions, we have from the equations 
(10) and (11) a differential equation in F as 



da 2 a"* 



Prom conformal transformation we have a = which implies ^ 



da _ (ln(/))m • 



m(f> 



(17) 



-. Now 



using equation (13) and (16), we have a transformed form of the equation (17) given by 



dY 2a; + 3 , ^ - 



4 ) 35o / a + 1 ^ ^ 
o{—r — — ln0 + Cln( 
3 [ 4 y n[a + 1) ^ 



-n(a+l) 



1 

a + l 



V'o"'(lnV')'"-i 2A'; ,^ ^ ,^ 

^ 7 + —{ln(p) r^-^ij,- 

m(p mcp 



(18) 

Using all the substitutions and the equation (9), (10) we get the expression for the potential 
V (after simphfication) as 



y = 



6-^^m\\n<t>f)Y + 6k^^'' 



00 



2(ln0)^V'o^ 1^(1 - " + \ ln0 + C(ln0)-"(--+^) 
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4 \y n{a + l) 
where Y is the solution of the equation (18). 



(19) 



The analytical function / in terms of is given by 



'In. 



(20) 



We now graphically represent V against (f) and / against 0. Prom the figure 1 we see 
that considering BD scalar field as Chameleon, the potential V{(f)) is sharply increasing 
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Fig.l Fig. 2 

Fig. 1 shows the variation of V against (p for flat (red curve), closed (blue curve), open (yellow 
curve) universe respectively and fig. 2 shows the variation of / against (p for 
a = 0.5, Bo = 4,C = l,m = 4,n = 2,V'o = 1,^^ = -1/3. 



with respect to the field and from figure 2 we see that the analytical function /(</>) is 
also increasing as field increases for some particular values of the parameters i.e., for 
a = 0.5, Bo = A,C = l,m = A,n = 2,ipo = l,uj = -1/3. 

III. DISCUSSIONS 

Here we have considered initially general scalar tensor Brans-Dicke theory of gravity having 
BD scalar field as Chameleon and then using conformal transformation we rewrite the field 
equations in Einstein's frame. Now considering the universe filled with variable modified 
Chaplygin gas and using some particular assumptions, we solve the density parameter for the 
radiation era. Now we graphically analyzed the potential of BD-Chameleon scalar field with 
respect to the field. We see from figure 1 that considering BD scalar field as a Chameleon and 
in presence of VMCG the potential is increasing with the increase of the field. We also see 
from the figure 2 that the analytical function of the field which is coupled with the matter 
Lagrangian in the action is also increasing with the increase of the BD-Chameleon scalar field. 
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